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We extend the model of a 2d solid to include a line of defects. Neighboring 
atoms on the defect line are connected by springs of different strength and different 
cohesive energy with respect to the rest of the system. Using the Migdal-Kadanoff 
renormalization group we show that the elastic energy is an irrelevant field at the 
bulk critical point. For zero elastic energy this model reduces to the Potts model. 
By using Monte Carlo simulations of the 3- and 4-state Potts model on a square 
lattice with a line of defects, we confirm the renormalization-group prediction that 
for a defect interaction larger than the bulk interaction the order parameter of the 
defect line changes discontinuously while the defect energy varies continuously as a 
function of temperature at the bulk critical temperature. 

PACS numbers: 05.10.Ln,05.10.Cc,62.20.-x 



I. INTRODUCTION 

Thermodynamics of solids with defects is a topic of current interest.-^ In this paper we 
expand an equilibrium statistical mechanics model^ of a solid to include extended defects. 
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Previously we usee™ the realistic anharmonic energy versus atomic distance developed by 
Ferrante and collaborators^ in a mean-field type computation. We then evaluated^ the role 
of thermal fluctuations by using renormalization-group and Monte-Carlo simulations. 

In our model described in Section II, the solid is constituted of harmonic springs. If the 
energy of such a spring is larger than a threshold, the spring is likely to fail.- In the limit of 
zero elastic energy the model reduces to the Potts model,- which has been used to describe 
correlated and uncorrelated percolation processes.-^ 1 ^^ 

In Section III we present phase diagrams based on the renormalization-group Migdal- 
Kadanoff scheme. This scheme is of course an approximation for models on regular lattices 
but as demonstrated by Berker— it is exact for hierarchical lattices. This latter feature 
and its simplicity make the Migdal-Kadanoff scheme quite popular.— il^iMiiZiM We analyze 
the stability of fixed points and find that the elastic energy is an irrelevant field in the 
renormalization-group sense. For this reason, in Section IV we report Monte-Carlo simu- 
lations of the model with zero elastic energy, i. e. Potts model on a square lattice with a 
defect line. Our goal is to verify the renormalization group prediction^ that for q = 3,4, on 
the defect line, at the bulk critical temperature, the order parameter jumps discontinuously 
from zero (high temperature phase) to a nonzero value (low temperature phase)while the 
energy varies continuously with temperature. This is interesting since for those q values 
the bulk transition is continuous. It was argued 19 that this unusual Id defect transition is 
due to the infinite range correlations at the bulk critical point. It resembles the Thouless 
transitions^ in a one-dimensional system with inverse squared distance decaying interac- 
tions. Our concluding remarks are found in Section V. 



The solid is made of "springs" some of which are live and some are failed upon ther- 
mal excitations. All processes are assumed to be reversible, unlike the work of Beale and 
Srolovitz22 where springs fail irreversibly. The energy of a "spring" < i, j > is given by the 
Hooke law: 



where r*j is the displacement vector from the equilibrium position of atom i, Ec is the 
cohesive energy, k is elastic constant and a is the equilibrium lattice spacing. If the energy 



II. MODEL 




(1) 



of the spring is larger than the threshold energy E the "spring" is more likely to fail than 
to be alive, p is the probability that the "spring" is alive and l — p the probability that the 
"spring" breaks. We assume its dependence on energy to be given by the Boltzman weight: 



_JL_ = = we -§ in-rl? (2) 

l—p 

E C + E 

where K = k/kgT and w = e k B T . 

For the extended line of defects the elastic, cohesive and threshold energies may take 
values different from the rest of the system. Hence while in the bulk the parameters are K 
and w, on the line of defects they are and Wd- 

We allow for correlations between failing events by using the Potts number of states q, 
which plays the role of a fugacity controlling the number of clusters. For q = 1 we have 
random percolation as springs fail independently. The partition function is a sum over all 
possible configurations of "live" springs: 

con jig 

C is the number of clusters, including single site clusters, and B is number of live "springs". 
The restricted partition function associated with the elastic energy for a given configuration 
of bonds (live " springs" ) is 

v config rr — elastic i a\ 

Elastic = Trve 3 ( 4 ) 

£ f^-f;) 2 (5) 
Kb± <i,j> A 

In Eq. © the sum is over all live " springs" . 

By using the Kasteleyn-Fortuin expansion^ for Potts model we can rewrite the partition 
function as 

Z = TrJTrre'^ (6) 

The Hamiltonian is 



B <i,j> ~ <i,j>defect ~ 

(7) 

where <jj is a Potts spin taking q values. This mapping is a Gaussian approximation valid 
when, on the right hand side of Eq. [71 the elastic energy is small compared to the first 
energy contribution. In our Monte Carlo simulations, we use the Hamiltonian in Eq. (7) for 
integer values of q. The coupling constants J\ and J2 are related to the original parameters, 
w and K, as follows: 



Ji = ln(l + w) (8) 

J ld = ln(l + w d ) (10) 

hd = Kd—^—r (11) 
Wd + 1 



III. RENORMALIZATION GROUP 

The Migdal-Kadanoff recursion equations^^ 5 - for two dimensions are obtained by assum- 
ing that each atom coordinate varies in the interval (-1/2, 1/2), where the equilibrium lattice 
constant is 1, and also using the Gaussian approximation (small elastic energy): 



w' = [l + U(w,K,q)] 2 - 1 (12) 

K'w' = K[l + U(w, K,q)\ U(w, K,q) (13) 

w' d = [1 + U(w, K, q)} [1 + U(w d , K d , q)} - 1 (14) 

K' d w' d = ^K d [l + U{w, K, q)]U(w d , K d , q) + ±#[1 + U(w d , K d , q)]U(w, K, q) (15) 



where 



w 2 erf( x fK/4:) 

U(w, K, q) = . W . (16) 

q x /K/n + V8w.erf( x /K/8) 



The above recursion equations represent the Gaussian approximation of the exact solu- 
tions for hierarchical lattices. The renormalization group flows in the bulk parameter space 
(w, K) are governed by the following fixed points at K = (pure Potts model): 

i. w = (non-percolating live "springs"), 

ii. w = oo (percolating network of live "springs"), 

iii. w = w c (Potts critical point). 

A stability analysis at the bulk Potts critical point, (K = 0, w = w c ) yields the two 
eigenvalues: 

i. the thermal eigenvalue: Ai (for the direction along the K = axis) is always larger 
than 1, meaning the w — w c is a relevant field; 

ii. The other eigenvalue associated with the flow along the w = w c line away from the 
pure model (K =0) is A 2 < 1 for all q. This means that there is a line of points in the 
(w, K) flowing into, and thus is in the same universality class as, the pure Potts critical 
point (w c , 0). 

The bulk phase diagram (FigJTJ), for any given q, in the (w, K) plane shows two phases: 

I. solid with mostly live springs, 

II. "crumbling" solid with failed springs, separated by a critical line in the universality 
class of the g-state Potts model. 

The defect fixed points are obtained from Eqs. (fT3|) and ffTol) after setting the bulk fields 
w and K at their fixed point values (bulk material is critical): w = w c , K = 0. All the fixed 
points are obtained for Kd = 0. We find fixed points at Wd = oo and at Wd = w c , for all 
q. For q > q% = 3 a third fixed point emerges from the one at Wd = oo at Wd = w*. At 
1 = Q2 — 6.8 (when exponent a = 0) the two finite Wd fixed points exchange position. 

The stability of the fixed points for perturbations of the defect fields (w^, K d ) is deter- 
mined by 

i. The eigenvalue Awd governing flow along Kd = axis is larger than unity (i.e. relevant 
perturbation) at the fixed point Kd = 0, Wd = max(w c , w*) and is less than unity (i.e. 
irrelevant perturbation) at the fixed point Kd = 0, Wd = min(w c ,w*) (Fig. EJ); 

ii. the eigenvalue AKd < 1 for all fixed points (Fig. [3]). 

The defect phase diagram (Fig. Hj) represented in the plane (w^, K d ) is obtained for the 
bulk fields at their critical values: w = w c , K = 0. 

The plane (wd, Kd) is divided in two regions by a critical line that flows into the fixed 
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FIG. 1: Bulk phase diagram for q = 10 and d 



2, in the plane (w,K). 



point: Wd = max(w c ,w*) , K d = 0. To the right of it (large wa), the live springs on the 
defect line are percolating, i.e. the Potts order parameter is nonzero on the defect line even 
though in the bulk it is zero. All points in that region flow to the fixed point Wd = oo, 
Kd = 0. The phase to the left of the critical line (small Wd) is governed by the fixed 



characterized by infinite correlation length or algebraic decay of correlations along the defect 
line and zero Potts order parameter. 

Since the elastic energy is irrelevant RG field, we concentrate next on the K = Kd = 
case. The RG analysis^ of this problem (two-dimensional Potts model with a defect line) 
predicts for 2 < q < 4 a hybrid phase transition on the defect line: discontinuous order 
parameter and continuous energy. 



point at Wd = min(w c ,w*), Kd = 0. This is a Berezinskii-Kosterlitz-Thouless^ 1 ^'^ phase 
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FIG. 2: Renormalization group eigenvalue Awd vs q: solid line at defect critical fixed point; dotted- 
dashed line at the non-percolating critical fixed point; dashed line at the percolating critical fixed 
point. Two vertical lines are q\ = 3 and qi = 6.8 (see text for comments). 

IV. MONTE CARLO SIMULATIONS 

In this section since we set the elastic couplings to zero in Eq. (J7|), we perform MC 
simulation of the Hamiltonian 



where J is the interaction parameter between the bulk nearest-neighbors (NN) and Jd is 
that between two NN on the defect line (J is equal to J\ of Eq. [HI and Jd is J\d of Eq. [TOl . 

We consider a square lattice of size N x x N y where N x = 80, 120, 160,200, 1920 and 
N y = 40, 60, 80. Each lattice site is occupied by a q— state Potts spin. We place the defect 
line at the middle of the y side of the sample, i.e. y = N y /2. The length of the defect line is 
thus N x . The reason why we use a very long length for the defect line stems from the fact 
that the determination of a transition in a linear chain needs large enough sizes to avoid 
statistical fluctuations. We use periodic boundary conditions. 

Our purpose here is to test the following RG prediction of the previous section: at the 
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FIG. 3: Renormalization group eigenvalue Al-Q vs q: solid line at defect critical fixed point; dashed- 
dotted line at the non-percolating critical fixed point; dashed line at the percolating critical fixed 
point. Since these eigenvalues are less than unity, the elastic energy is an irrelevant field at the 
three fixed points. Two vertical lines are q\ = 3 and q2 = 6.8 (see text for comments). 

bulk transition temperature T c , the defect line undergoes a phase transition where its order 
parameter is discontinuous but its energy is continuous. The bulk transition temperature is 
the temperature at which the phase transition of the system without the defect line takes 
place. 

Let us consider the case where q — 4. The critical temperature is given by the exact 
formula {k B T c / J)' 1 = ln(l + y/q). With q = 4, one has T c ~ 0.910239 in units of J/k B . 
Note that this value of T c corresponds to the thermodynamic limit, i.e. infinite system 
size. In MC simulation, we work at finite sizes, so for each size we have to determine the 
"pseudo" transition which corresponds in general to the maximum of the specific heat or of 
the susceptibility. The maxima of these quantities need not to be at the same temperature. 
Only at the infinite size, they should coincide. The theory of finite-size scaling permits to 
deduce properties of a system at its thermodynamic limit. We have used in this work a size 
large enough to reproduce the bulk transition temperature up to the fourth decimal. 




FIG. 4: Phase diagram at w = w c , K = in plane (w^,K^). The two critical phases (infinite 
correlation length) are separated by a defect critical line. In the small Wd phase there is zero 
probability for percolation on the defect line, while in the large Wd phase there is a finite probability 
for percolation on the defect line. 

In order to determine the nature of the phase transition of the defect line, we shall use 
the histogram technique^ which is known to allow us to distinguish with accuracy the order 
of the phase transition. 

The simulation is carried out as follows. We fix J = 1 hereafter. For each value of 
Jdi using first the standard Metropolis MC method^ we equilibrate the system of a given 
size N x x N y at a given temperature T during 10 6 Monte Carlo sweeps (MCS) per spin 
before averaging physical quantities over the next 2 x 10 6 MCS. We determine the transition 
temperature at the given size N x x N y by examining the calculated physical quantities such 
as the internal energy per spin E, the specific heat C v per spin, the Potts order parameter 
Q and the susceptibility per spin \- F° r the bulk g-state Potts model, Q is defined as 

qmax(Q 1 ,Q 2 ,...,Q q ) - 1 , , 

Q- — (18) 

where Qi = N n ^ N {i = 1, q), rii being the number of sites having qi. For the defect line, 




0.82 0.84 0.86 0.8 



- v . 



0.B4 0.86 0.8. 



0.92 0.94 0.9i 



FIG. 5: Case without defects: Energy per spin E (upper curve) and specific heat per spin C v 
(lower curve, line is a guide to the eye) vs temperature T for q = 4, with N x = 1920, N y = 60 and 
Jd = J = I- 

the order parameter Qd is similarly defined on the defect line, namely 



Qa 



gmax(Q / 1 ,g / 2 , ...,Q' q ) - 1 
9-1 



(19) 



where Q- = ^- [i = 1, g), being the number of sites on the defect line having g;. 

Let us show first in Fig. [5] the energy and the specific heat of the case without defects 
where the bulk critical temperature is T c ~ 0.9103 for q = 4 with the size used here (the 
exact value of T c is 0.910239). The energy of the defect line is shown in Fig. 

The order parameters the bulk and of the defect line are shown in Fig. [7] for q — 4, J — 1 
and Jd = 1.6 and 2. 

These figures show that at the "bulk" phase transition temperature T c = 0.9103, while 
the defect energy is continuous, the order parameter of the defect line undergoes a vertical 
fall, indicating a discontinuity predicted by the RG analysis shown in the previous section. 
For large values of Jd, for instance Jd = 2, the defect line takes a very long time to become 
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FIG. 6: Energy per spin Ed of the defect line vs temperature T for q = 4, Jd = 1.6 with N x = 1920, 
N y = 60 and J = 1. Note that Ed is continuous at T c . See text for comments. 
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FIG. 7: Order parameters Qd of the defect line for Jd = 2 (upper curve with black circles), 
Jd = 1-6 (middle curve with void circles) and bulk order parameter Q (lower curve with crosses) 
vs temperature T for q = 4 with N x = 1920, N y = 60 and J = 1. Note that Qd is discontinuous at 
T c . See text for comments. 



disordered for T slightly larger than T c . Several millions MC sweeps are necessary. 

In order to check the behaviors of E d and Qd at T c , we have calculated the defect energy 
histogram not only at T c but also in the temperature region around T c . As it turned out, we 
observe only a gaussian distribution of Ed (see Fig. [8] confirming the absence of discontinuity 
of Ed. Note that if the energy is discontinuous, its histogram should show a double-peak 
structure, not a gaussian one. 

We have also established a histogram for Qd in the following manner. We divided the 
interval between and 1 into N x intervals. At each MC sweep, we added 1 in the interval 
corresponding to the value of Qd- In doing so for 2 millions MC sweeps, we obtained 
a histogram for Qd which is shown in Fig. [9] at T c . As seen, we have a double-peak 
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FIG. 8: Energy histogram of the defect line taken at T c , for q = 4, Jd = 1-6 with N x = 1920, 
N y = 60 and J = 1. See text for comments. 




FIG. 9: Histogram of the order parameter Qd taken at T c , for q = 4, Jd = 1.6 with N x = 1920, 
N y = 60 and J = 1. Note the double-peak structure indicating that Qd is discontinuous at T c . See 
text for comments. 

distribution of Qd, indicating that during the time the defect line can have both the ordered 
and disordered phases. This is a strong signature of the discontinuity of Qd- 

Note that the size effects for N y = 40, 60, 80 and 100 are not significant and are included 
in the error estimation. Simulations have been carried out also for other values of Jd between 
1.2 and 4. The results show the same aspects as those shown above with Jd = 1.6. 

Now, let us examine the case where q = 3. The results are very similar to the case q = 4 
shown above. So the conjecture of the RG analysis where the transition of the defect line 
is of first order for J q > J at T c for q > 3 is verified here. Figure \IU\ shows the defect order 
parameter Qd versus T for q = 3 with several values of Jd- Note that T c = 0.9949 for q = 3. 

Let us show in Fig. [TT]the value of Qd taken at T c as a function of Jd- As seen, the gap 
is increased with increasing Jd- 



FIG. 10: Order parameters Qd of the defect line for Jd = 2 (upper curve with black circles), 
Jd = 1.4 (middle curve with void circles) and bulk order parameter Q (lower curve with black 
squares) vs temperature T for q = 3 with N x = 1920, N y = 60 and J = 1. Note that Qd is 
discontinuous at T c . See text for comments. 
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FIG. 11: Value of the order parameter Qd taken at T c versus Jd for q = 4 (upper curve) and q = 3 
(lower curve) with N x = 1920, N y = 60 and J = 1. Lines are guides to the eye. 

At this stage, it is worth to mention that due to the usual finite size effect, the order 
parameters Q and Qd do not vanish above T c : a finite tail exists and decreases with in- 
creasing lattice size. For Q, the transition point is, in simulations, taken at the change of 
curvature of Q, i.e. at the maximum of the corresponding susceptibility. As for Q d , due to 
its discontinuity at T c , the values shown in Fig. [TT]are the upper one at T c . For clarity, the 
error bars are not shown there, but it is on the second digit, for instance Qd = 0.880 ± 0.020 
for Jd = 1.6. 

We show now the time dependence of the order parameter Qd and the energy. Figures [12] 
and H3] show these quantities for Jd = 2 at T c = 0.9949 and at a temperature slightly above 
T c , namely T = 1. Two procedures have been used: i) heating, i.e. using the ordered phase 
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FIG. 12: Time dependence of the energy par spin of the defect line for q = 3 with Jd = 2 at 
T c = 0.9949 and at T = 1. From the bottom: heating to T c , cooling to T c , heating to T = 1, 
cooling to T = 1. See text for comments. 

(F) as initial spin state ii) cooling, i.e. using a paramagnetic state (P) as initial condition. 
We discuss first the energy case. At T c , heating from the F state (bottom curve in Fig. 
IT~2T) and cooling from the P state (second curve from the bottom in Fig. fl"2l) give the same 
energy only after two millions MC steps/spin (see the first million MC steps in Fig. fl~2l) . 
At T = 1, one needs almost the same MC time to get the same energy for the heating and 
cooling (curves three and four from the bottom in Fig. [T2l) . 
The time dependence of Qd is interesting (Fig. [lH]) : 

i) At T c , Q d with F condition stays stable (first curve from top in Fig. [T3]) and reaches 
the stationary value 0.88 at two millions MC steps/spin. This means that the defect line is 
ordered at T c while the bulk spins are disordered. However, when cooled from the P state 
(third curve from the top), Qd stays very small (disordered state). We conclude that at T c 
there are two possible values of Q d for the same energy. This explains the gap of Qd at T c 
shown earlier. 

ii) At T = 1, slightly above T c , Q d takes a long time (~ 500000 MC step/spin) to become 
disordered in the heating procedure (second curve from the top) while it is disordered all 
the way in the cooling procedure. 

Let us discuss more about the discontinuity of the order parameter Qd at T c . In a standard 
first-order transition, the ordered and disordered phases coexist at T c yielding a discontinuity 
in both energy and order parameter. In the case studied here, only Qd is discontinuous, while 
Ed is not. Thus, there is no double-peak energy distribution. To answer the question how Qd 
can take two values at T c , we have examined the snapshots of the defect line taken during 
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FIG. 13: Time dependence of the order parameters Qd of the defect line for q = 3 with Jd = 2 at 
T c = 0.9949 and at T = 1. From the top: heating to T c , heating to T = 1, cooling to T c , cooling 
to T = 1. See text for comments. 

the simulation time. Interesting enough, we observed configurations schematically of the 
types: 

i) 111111. .1111111. .111111.. where dots indicate Potts values other than 1 (same kind of 
ordered segments separated by disordered portions) 

ii) 111111. .3333333. .222222.. (random ordered segments separated by disordered por- 
tions) 

It is obvious that these configurations give the same energy if the number of disordered 
portions is the same, while Qd is different: the first configuration yields a large value of Qd 
and the second a zero one. This observation explains why at T c the energy is continuous 
but the order parameter is not. 

The phase transition observed here for the defect line is very interesting in several aspects: 

i) the nature of the transition is novel in the sense that only the defect order parameter 
is discontinuous, not the defect energy, nor the bulk order parameter and bulk energy; 

ii) note that long-range interaction in one dimension can cause a first-order 
transition,— iM^l while systems with short-range interaction do not show such a phase tran- 
sition. The present defect line with NN interaction shows thus an exception. We believe 
that its immersion in a disordered neighboring lines at T c plays a key role in provoking such 
a pseudo-discontinuous transition. 



V. CONCLUSIONS 



In this paper we have studied a model of a two-dimensional solid with an extended defect 
line using the Migdal-Kadanoff renormalization group. Since the elastic energy turns out to 
be an irrelevant field, we then studied this model for zero elastic energy, i.e. 2d Potts model 
with a defect line, using Monte Carlo simulations. The renormalization-group analysis of the 
Potts model on a 2d hierarchical lattice with a defect line, suggests an interesting behavior 
at the bulk transition T c if the interaction on the defect line J d is larger than the bulk 
interaction J: the order parameter of the defect line should undergo a discontinuity at T c 
while the energy is continuous. Our Monte Carlo simulations of the 3- and 4-state Potts 
model on the square lattice with a defect line confirm the RG prediction. 

To conclude, we would like to emphasize two points of interest. First, the existence of a 
phase transition in a one dimensional system (defect line) is rare. It is induced by the bulk 
infinite correlations at criticality. Second, the nature of the phase transition is unusual in the 
sense that the order parameter is discontinuous while the energy is not. This is somewhat 
similar to the Thouless transition in Id models with long-range interactions. 

M.K. wishes to thank the University of Cergy-Pontoise for hospitality while this work 
was carried out. 
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